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In this talk results of the study of the —)■ 7i~e^u decay channel, 

recorded by the BABAR detector at the c.m. energy close to 10.6 GeV, 
are reported. The branching fraction of this channel is measured relative 
to the —)■ K~n^ decay. The hadronic form factor, function of 

the four momentum transfer squared between the D and the tt mesons, 
is compared to various theoretical predictions, and the normalization 

= 0) is extracted from a fit to data. Results are compared with 
Lattice QCD calculations. A new multipole model is applied which makes 
use of present information of resonant states contributing to the form 
factor. With the understanding of the £)(g^) form factor, and provided 
the relation between the —)■ 7r~e^u and B^ —)■ 7i~e~^u decay widths 

at the same pion energy, the CKM matrix element V^b is determined and 
compared to recent measurements. This method of extracting I4b will 
become competitive with new Lattice QCD calculations of the ratio of 
form factors. 
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1 Motivation 


The differential decay width of the —)■ TT~e'^v decay channej^as function of the 
four momentum transfer squared between the D and the vr mesons, can be expressed 
in terms of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, Vcd, and, ne¬ 
glecting the electron mass, a unique form factor, • This form factor describes 

all the non-perturbative QCD effects in the D —)■ vr transition: 
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p* being the pion momentum in the rest frame. The form factor can be represented 
as an inhnitive sum of poles, corresponding to resonant states, D*, with = 1~ 
which couple to Dtt. 
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For the —>■ TT~e^v decay channel one has the advantage that part of fd. ^iq^) is 

known, since contributions from the leading state D* and the first radially excited 
state D* are known and can be used to constrain the form factor BE]. Another 
interest comes from the fact that the —)■ 7r“e+z/ and —)■ 7r“e+z/ decay channel 

can be related at the same pion energy, allowing, if the form factors are known, the 
extraction of the CKM matrix element Vub- 


2 Analysis method 

This analysis [I] is based on similar techniques as other charm semileptonic decays 
at BABAR [3l IU |5]. The —)■ 7r~e~^iy decay channel has the difficulty that it is 

Cabibbo-suppressed, with a small branching fraction, and suffers in addition from 
large background from pions. Using 347.2 fb“^ of e’''e“ —)■ cc data recorded by the 
BABAR detector at the T(4S') energy, the decay D*~^ —)■ D^7r~^ with —>■ 7r~e~^i' 

is reconstructed using a partial reconstruction technique. The two pions and the 
positron are reconstructed in the same event hemisphere, requiring a tight particle 
identification for signal pions and vetoing kaons. The four momentum is obtained 
from the reconstructed particles and the missing energy from the information of the 
rest of the event. Constraints on the D and D* masses are applied in a kinematic fit 
to obtain the distribution. This method is validated using hadronic —>■ 

data. The main issue of the analysis concerns the background suppression. Using 
Fisher discriminant variables against BB and charm backgrounds, the S/B rate is 

* Charge conjugated are implicit in this document 
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about 1.2, with a signal efficiency around 1.8%. To further control the several sources 
of background events, the mass difference between the D* and the is used. The 
signal region is dehned as 6{m) < 0.155 GeV^. Two additional 6{m) windows are used 
to evaluate the background contributions. Using the event missing energy and the 
pion momentum information, the rates for the different types of background events 
are constrained. In this way, the main source of systematic uncertainty in the analysis, 
due to the background control, is assessed using data. The 6{m) distribution is shown 
in Fig. left. The distribution, = (p^o — = {Pe+ + Pu)'^ is measured then 
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Figure 1: Left: mass difference 6{m) = m{D^7r^) — m{D^). Open and full circles 
correspond to kinematic hts applying the and D* constraints pQ. The distribution 
for MC-simulated signal and the different background distributions are superimposed. 
Right: the measured g^ distribution (data points) for events selected in the 6{m) 
signal region is compared to the sum of the estimated backgrounds and the htted 
signal component. 


in 10 bins and a £t to data is applied using several parameterizations of the form 
factor. In Fig. [fright the g^ data distribution is shown. The MC-simulated events, 
with corrected background components is also shown. It is verihed that the angular 
distribution is well reproduced for the signal and background components. 

—)■ 7r~e~^i> events are normalized to the number of K~'k^ decay events, 

trying to have a selection as similar as possible for the semileptonic and hadronic 
channels. The ratio Rr, = is measured to be Rd = 0.0702 ± 0.0017 ± 

0.0023, where the hrst uncertainty is statistical and the second one is systematic. 
Using the world average for the B{D^ —)■ K~n^) [6j, the branching fraction of the 
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—)■ 7r“e+i/ decay channel is measnred to be —)■ 7r“e+n) = (2.770 ± 0.068 ± 

0.092 ± 0.037) X 10“^, where the third uncertainty comes from the uncertainty in the 
normalization channel. 


3 Form factor interpretation 

Having measured the number of —)■ 'K~e'^v events as function of (Fig. [fright), 
several parameterizations can be tried to describe and £t the form factor. One of the 
most extensively used is called the 2 ;-expansion, a model-independent parametrization 
based on general properties of QCD [7]. In terms of the variable 2 ;, dehned as 


■\/t+ — t + — to 
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where t = q^, to = t+(l — — t-/t+) and t± = {rrijjo ±mjr+)^, the form factor, takes 

the form: 

1 00 

. A E«fc(^o) (4) 

01 fc=0 

where P{t) = 1 and $(f,fo) is an arbitrary analytical function. A standar choice for 
<F(t, to) is used by the different experiments pQ. The htted parameters are commonly 
dehned as = a^/ao for k = 1,2, and the overall normalization of the expansion 
is Vcd X /+ £)(0). Results of the £t to BABAR data are: ri = —1.31 ± 0.70 ± 0.43, 
ra = -4.2±4.0±1.9, and KdXflniO) = 0.1374±0.0038±0.0022±0.0009. The last 
uncertainty in the normalization comes from external inputs [T] . The main disadvan¬ 
tage of this parametrization is that the parameters have no physical meaning and 
cannot provide an interpretation of the form factor. It is difficult to constrain the 
contribution from the D*~^ pole in this parameterization because it requires extrap¬ 
olation beyond the physical region. Fig. [^compares the BABAR result with several 
measurements at different experiments and with Lattice QCD calculations. Assuming 
unitarity of the CKM matrix, with \Vcd\ = jlAsI = 0.2252 ± 0.0009 |S], the BABAR 
result leads to /+£>(0) = 0.610 ± 0.017 ± 0.010 ± 0.005. If, instead, one uses the 
Lattice QCD results on the form factor, /+ £)(0) = 0.666 ± 0.029 |9], the Vcd matrix 
element results in Vcd = 0.206 ± 0.007 ± 0.009. 

However, one can go further in the understanding of the —)■ 7r~e~^u form factor: 
based on HH and [12], a three-poles ansatz has been developed P]. The form factor 
can be expressed as an inhnitive sum of = 1~ states (Eq. |^. The residue which 
defines the contribution of these states {D* resonances) can be, in turn, expressed in 
terms of the meson decay constants fn* and coupling to the Dn hnal state, gD*Dn- 

= ^mD;fD;gD^- (5) 
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Figure 2: Left: form factor fit using the z-expansion for BABAR [T], CLEO-c [lO] and 
the HFAG average [6]. Bands represent the £t uncertainty. Right: Comparison of 
/+d( 0) obtained by different experiments and by Lattice-QCD calculations. 


For the —)■ 7r~e~^i> decay channel, the decay constants for the two hrst states, 

the leading D* meson and the hrst radially excited state, D* , have been computed 
by Lattice QCD |2]. The couplings can be obtained from the experimental 

results at BABAR of the measured widths for these resonances [131 [H]. fhis way the 
contribution of these states to the form factor is well dehned and constrains 
In Figure [^left the contributions of the D* and the D*' are given, together with the 
form factor £t to BABAR data using the 2 ;-expansion formalism. This plot reveals the 
fact that the form factor cannot be explained by only these two contributions, and 
then additional poles are needed to £11 the gap between the measured data and the 
D* + D* contribution (gray curve). Using the constraints given by the D* and the 
D* poles, one can consider an ehective third pole contributing to the form factor [2] . 
The superconvergence condition, — 0 , is applied, from the behaviour 

of the form factor at very large values of In Figure [fright, the £t to BABAR 

data using this description of the form factor is presented. Once should note that 
data is well described by this ansatz. The e£ective third pole mass is htted and 
results in the value = (3.6 ±0.3) GeV , which is larger than the predicted third 

e// 

= 1~ state by quarks models (~ 3 GeV), as it is expected since it is considered 
as an ehective pole. A unique contribution from this predicted state is excluded by 
data. 
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Figure 3: Left: Measured /+£,(g^) at BABAR and expected contributions of the 
D* and D*' poles to the form factor. Bands represent the uncertainties. The gap 
between data points and the gray curve reveals the need of additional hadronic states 
to describe the form factor. Right: BABAR data and htting curves using the z- 
expansion (solid blue) and the three-pole ansatz (dashed green). 


4 Vub extraction 

Having measured dT, one can extract the CKM matrix element V^b from 
the relation between the —)■ and B^ —)■ 7r“e+z/ decay channels, valid for a 

common range in the energy of the ejected pion in the rest frame of the heavy-light 
meson, which is less than about ~ 1 GeV. Instead of E*, one can use the Lorentz 
invariant variable wd(b) = vd{b) -v-n-, where vd{b) = Pd{b)/^d{b) and = p^/rriT, are 
the four-velocities of the D{B) and tt mesons, respectively. In terms of this quantity: 
nil — 2 mD{B)mT,WD(B), and then at Wd = Wb- 


dT^/dwB rrtB fVu, 


dT^/dwD mo 


(Vuby 

fl,BiwB) ^ 
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The Vub element can be obtained from Eq.(|^ if the ratio between the /+_s(ic_b) 
and /+£)(tCD) form factors is known. This ratio can be obtained from Lattice QCD 
calculations, or using a phenomenological model. The experimental common range 
for the two decays on wb,d is between 1 and 6.7. This corresponds to from 0 to 
2.975 GeV^ for D decays and from 18 to 26.4 GeV^ for B decays, leading only about 
17% of overlapping region. Nevertheless, a physics interpretation of the charm form 
factor allows to use it outside the D physical region and to determine the form factor 
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ratio. 

Two different approaches have been used to extract Vub, leading to systematic 
uncertainties of very different origin. BABAR data for —)■ 7r~e'^i> [T] and B^ —)■ 

7i~e~^h' decays [U] are used: 

- Vub from Lattice results of individual form factors: 

Considering the Lattice results for [Ml HZj and /+£)(g^) [M] one observes 

that the two form factors have a similar w dependence. For w > 4, the ratio of these 
form factors is 1.8 ± 0.2[^ Assuming this constant value, and extrapolating to the 
nonphysical region the measured —>■ 'n'~e~^u form factor, in terms of the three-poles 
model, the Vub matrix element is htted, as it is shown in Fig. |^left, giving the value: 

Vub = (3.65 ± 0.18 ± 0.40) x lO-^. 

The dominant contribution to the systematic uncertainty originates from the form 
factor ratio from Lattice. This result can be improved if Lattice QCD provides values 
for this ratio with better accuracy and for several values of g^. 

- Vub from the three-poles phenomenological model: 

Another approach consists in htting directly B^ —)■ 7r~e~^i> events using the three-pole 
phenomenological model for the B ^ tt form factor, since it has been proven to work 
well for the —)• 7r~e'^i> decay. For that one can use constraints from the residues 

for the hrst two poles, which correspond to the B* and B^ , taking the mass of the 
latter from [191, and htting the third pole with an effective mass. The result of this 
third effective pole is = (7.4 ± 0.4) GeV. It is expected that the ratio for the 

e// 

residue at the different poles are the same for D and B semileptonic decays nn, and 
a constraint is applied in the ht. The ht results in: 

Vub = (2.6 ± 0.2 ± 0.4) X 10-3, 

and it is shown in Figure [fright. The uncertainty on Vub is dominated in this case 
by the knowledge of the gB*(.')B-w couplings entering in the residue. This value could 
also be improved by new Lattice calculations. 

5 Conclusions 

The —)■ Ti~e'^v form factor and branching fraction is measured at BABAR . Re¬ 

sults are competitive and in agreement with CLEO-c, BELLE and preliminary results 
from BES III. Experimental results in this channel are at present more accurate than 
Lattice QCD calculations. A physics interpretation of the D —)■ vr form factor is de¬ 
veloped, in terms of a phenomenological three-poles model, using precise information 

^One should note that this value is consistent with the expectation at first order: = 

1 . 7 . 
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Figure 4: Left: Fit to the dB/dw differential decay rate for —)■ TT~e^v events 
measured by BABAR with an extrapolation of the —)■ 7r~e~^i> form factor mea¬ 
surement, assuming a constant ratio of form factors equals to 1.8 ± 0.2 from Lattice 
calculations. The Vub parameter is htted. Right: Fit to the —)■ 7i~e^h' events using 
the three-poles phenomenological model for the form factor, using constraints for the 
residues of B*,B* and htting a third pole with an effective mass. The contribution 
for each pole is shown. The two lines indicate the theoretical uncertainties. 

of the two hrst contributing poles. It is observed that these two poles cannot explain 
alone the form factor and an effective third pole is obtained. This description agrees 
well with data. The Vub matrix element is extracted using charm semileptonic de¬ 
cays, through two alternative approaches: assuming a constant form factor ratio from 
Lattice QCD, or considering the three-poles model, proved on —)■ decays, 

on B^ —)■ n~e~^i' events. Results are in agreement with recent Lattice calculations 
m and the measurement from LHCb EH. in particular for the hrst approach. This 
method of Vub extraction will become competitive if new Lattice QCD computations 
on the form factor ratio are available. 
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